Progressive suppression of spin relaxation in 2D channels of finite width 
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We have investigated spatio-temporal kinetics of electron spin polarization in semiconductor nar- 
row 2D strip and explored the ability to manipulate spin relaxation. Information about spin of the 
conduction electrons and mechanisms of spin rotation is incorporated into transport Monte Carlo 
simulation program. A model problem, involving linear-in-fc splitting of the conduction band, re- 
sponsible for the D'yakonov-Perel' mechanism of spin relaxation in the zinc-blende semiconductors 
and heterostructures, is solved numerically to yield the decay of spin polarization of an ensemble of 
electrons in the 2D channel of finite width. For very wide channels, a conventional 2D value of spin 
relaxation is obtained. With decreasing channel width the relaxation time soares rapidly by orders 
of magnitude. Surprisingly, the cross-over point between 2D and quasi- ID behavior is found to be 
at tens of electron mean-free paths. Thus, classically wide channels can effectively suppress electron 
spin relaxation. 

PACS numbers: 72.15.Lh; 76.60.Es 



I. INTRODUCTION 

SpintronicsJj'EI a nontrivial extention of the conven- 
tional electronics, adds functionality utilizing the carrier 
spin degree of freedom. Spin can potentially be used 
as a by far more capacious quantum information stor- 
age cell, be involved in the transfer of information, for 
elaborate schemes of information processing, both quan- 
tum and classical, and be integrated with electric-charge 
counterparts in combined designs. Electron or nuclear 
spin manipulation is believed to be the key component 
of the potential realizations of a quantum computer (for 
example, see Ref. ||). 

Several devices were proposed so far including spin 
hybrida and field effecip.transistors,Q tunneling structures 
with magnetic layers,l3 and spin-based memory.Q The 
most notable success was obtained with giant magnetore- 
sistance (GMR) effect devices that rely on the variation of 
electron scattering in a multilayer stack oLferromagnetic 
films separated by nonmagnetic materialsH'u Though the 
variation of the current through the structure is small, it 
is sufficient for detection, and the excellent sensitivity 
of the device to weak external magnetic fields (typically 
1% change of resistance per oersted) opened the door for 
massive data storage applications. Spin coherence is a 
pivotal prerequisite for the operability of the prospective 
spintronic devices. 

Concerning experimental realization of spintronic el- 
ements, most propositions rely on the injection of 
spin-polarized electrons from the ferromagnetic layer 
and suffer from the poor quality of the ferromagnet- 
semiconductor interfaces that produce a large number of 
surface states causing strong spin relaxation and reduc- 
ing the polarization of the injected electrons to only a few 
percent. Since the injected electron crosses the metal- 
semiconductor interface twice in spin- valve designs, i.e., 
at the source and drain terminals of the device, the im- 
portance of the interface problem multiplies. Further- 
more, spin relaxation in the active region of the device 



adds to the complications as well. 

In our study, we consider the possibilities for suppress- 
ing spin relaxation of the conduction electrons and eval- 
uate different approaches. It is extremely interesting to 
investigate what happens with the spin relaxation in a 
2D electron gas as the long strip of finite width is formed 
using electrostatic squeesing split-gate technique or by 
another method. These systems can be used to connect 
active elements in integrated chip designs and can even 
become a part of the active device as we continue to re- 
duce size of the element. Eventually we should reach the 
ID limit exhibiting no spin relaxation. This transition 
to the ID behavior was recently documented in the com- 
puter simulation of the Diatta and Das spin transistor 
reported by Bournel et al.ll3 that inspired our research. 
We concentrate our efforts on the definition of the cross- 
over regions from the point of view of the spin relaxation 
time and the description of its behavior in the broad re- 
gion around the cross-over points. 

The rest of this paper is organized as follows: We start 
with a brief account of spin relaxation mechanisms in 
semiconductors in order to identify the most relevant. We 
discuss its transformation to lower-dimensional systems, 
as well as different possibilities to suppress destructive 
spin relaxation (Sec. O). The model of a narrow pat- 
terned 2D electron gas channel is described in Sec. III. 
Results of Monte Carlo simulation of spin relaxation in 
the channel are given in the next section (Sec. IV), along 
with a comprehensive discussion of the relaxation regimes 
in terms of the channel width and a value of the spin 
splitting of electron subbands (Sec. |V|). A brief summary 
follows at the end. 



II. MECHANISMS OF SPIN RELAXATION 

There are several mechanisms-that can cause spin re- 
laxation of conduction electronsliil (see Ref. 11^ for an up- 
to-date informal review): 



i) The mechanism of D'yakonov and Perel' (DP)ll3 takes 
into consideration that spin sphtting of the conduction 
band in zinc-blende semiconductors at finite wave vectors 
is equivalent to an effective magnetic field that causes 
electron spin to process. For electron experiencing ran- 
dom multiple scattering events the orientation of this ef- 
fective field is random and that leads to the spin relax- 
ation. l— I 

ii) Bir-Aronov-Pikus (BAP)ll3 processes involve a simul- 
taneous flip of electron and hole spins due to electron-hole 
exchange coupling. 

Hi) Spin relaxation due to momentum relaxation is pos- 
sible |-d|irectly through spin-orbit coupling [Elliot- Yafet 
(EY)y process]. 

iv) Spin relaxation can take place as a result of hyper- 
fine interaction of electron spins with magnetic momenta 
of lattice nuclei, the hyperfine magnetic field being ran- 
domly changed due to the migration of electrons in the 
crystal. 

BAP processes require a substantial hole concentration 
that is not available in the unipolar doped structures. EY 
processes are suppresed in 2D environments — target of 
our prime interest. Thus, we concentrate our efforts on 
the DP mechanism as the most relevant one for the case 
considered. 



it is rather easy to design a structure with \kz\ > kp, 
the splitting for a typical electron will be larger and the 
rate of spin relaxation will be enhanced. Similar behav- 
ior can be naively expected with the further reduction 
of the dimensionality to the ID case. Let the axis x be 
along wire in what follows. Here again the main terms 
will contain spacial-confinement multiplyers ky, k^. The 
principal difference with the 2D case is that now all ro- 
tations are limited to a single axis direction and they 
commute with each other. Apart from the systematic 
rotation, spin polarization does not disappear with time. 
All particles, independently of the number and the se- 
quence of the scattering events, that reach the same final 
point B will have the same spin orientation. This state- 
ment is relaxed if one allows intersubband scattering in 
the one-dimensional system. This type of scattering be- 
comes progressively more important for wider and wider 
quantum wires with more and more subbands involved. 
Thus, we will recover, as one can predict, the 2D or 3D 
case in the limit of very thick quantum wires. 

Actually, there are two types of terms that appear 
in the effective-mass Haniiltonian for the 2D electron 
gas: the bulk-asymmetry-I^TI and structure-asymmetry- 
induced (Bychkov-Rashba)Jl30 They are of the same 
functional form 



A. DP mechanism and change of the dimensionality 

As a result of the relatively low zinc-blende crystal 
symmetry, the effective 2x2 electron Hamiltonian for 
the conduction electrons contains spin dependent terms 
that are cubic in the electron wave vector k 



H = r]'[axky{kl-kl)+ etc 



(1) 



The constant rj' refiects the strength of the spin splitting 
in the conduction band whose value is defined by the de- 
tails of the semiconductor band structure, at {i = x, y, z) 
are the Pauli matrices, other terms in Eq. (m should be 
obtained by the cyclic permutation of the indices. 

As the system dimensionality is changed from the 3D 
to 2D by, for example, the extremely strong spacial con- 
finement in the third direction — that can be achieved in 
the semiconductor heterostructurealj — modifications to 
the character of the spin relaxation occur as well. First of 
all, an average wave vector in the direction of the quan- 
tum confinement (axis z) is large, so the terms in the spin 
splitting involving k^ will dominate. This results in the 
orientation of the effective magnetic field in the plane of 
the 2D electron gas {x — y plane). Still the elementary 
rotations around random axes, all laying in one plane, do 
not commute with each other, so the electrons reaching 
the same final destination by different trajectories will 
have different spin orientations. For longer times, more 
and more distinguishable trajectories will become possi- 
ble and this will lead to a progressive reduction of the av- 
eraged spin polarization of the electron ensemble. Since 



H = rj (T ■ [k X z] = rj{axky 



<^y l^x ) 1 



(2) 



where £ is a unit vector in the z direction. We are 
not interested here in the origin of the spin splitting 
term in our model problem and simply assume its pres- 
ence with a wave- vector dependence given by Eq. (y). 
(see Refs. 20 2^ for the band-structure calculation of the 
linear-in-fc spin splitting in heterostructures simultane- 
ously treating bulk- and structure- induced asymmetry.) 
This form of the spin Hamiltonian is equivalent to the 
precession of the spin in the effective magnetic field 



H 



n. 



off, 



where Jlcff — ?7dp v x z. 



(3) 



Here the particle velocity v — hk/m*\ and an obvious 
substitution rj — > rjDp is done for convenience, r/pp is 
expressed in inversed length units. For a particle moving 
ballistically the distance l/rjop spin will rotate to the 
angle (/) = 1. We remind that the quantum-mechanical 
description of the evolution of the spin 1/2 is equivalent 
to the consideration of the classic momentum 5' with the 
equation of motion 



dS_ 
dt 



= ricff X s. 



The reciprocal effect of electron spin on the orbital mo- 
tion through spin-orbit coupling can often be ignored due 
to the large electron kinetic energy in comparison to the 
typical spin splittings and strong change of the momen- 
tum in scattering events. 



B. Possiblities to influence spin relaxation time 

In addition to understanding of the reasons govern- 
ing spin depolarazation of carriers, we wish to consider 
and assess the possibihties to actively influence these de- 
structive processes in order to improve parameters and 
gain new functionahty of the future spintronic devices. 
FoUowing are the potential approaches for manipulating 
spin relaxation times: 

i) A rather simple observation follows directly from the 
essence of the regime of motional narrowing (small e\£^ 
mentary spin rotations during ballistic electron flights). EJ 
Since Tg ^ Tp(51^g), reduction of momentum relaxation 
time, Tp, leads to the suppression of spin relaxation. On 
the other hand, this will lead to increasing of broadenings 
as well as decoherence, and can worsen device parame- 
ters. 

ii) Since the bulk-asymmetry- and structure-asymmetry- 
induced spin splittings are additive with the same k- 
dependence, it is possible to tune combined spin split- 
tings in the conduction band to a desired value through 
manipulation of the external electric field. 
Hi) Additional spin splitting, which is independent of the 
electron wave vector will fix the precession axis. An ev- 
ident possibility here is the Zeeman effect. The time of 
spin relaxation scaJes in the presence of the external mag- 
netic field, B, asE3 
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1 



TsiB) r5(0)l + (r!LTp)2' 

where TifiL — gi^BB is a Zeeman splitting of electron spin 
sublevels. This equation suggests that for fii^Tp — 1 spin 
relaxation time will double. 

iv) There is a possibility of controlling the spin relaxation 
of the conduction electrons by doping. The flrst-jealiza- 
tion was reported in (S-doped heterostructures£3 and a 
several-orders-of-magnitude longer spin metcupry has re- 
cently been observed in n-doped structures .cfl 
v) When the channel width, L, is comparable to the mag- 
nitude of the electron mean free path, Lp, from the clas- 
sical point of view, the sequential alteration of one of the 
wave vector components should e|fiEectively reduce spin 
relaxation (reflective boundaries) JlJ Scattering on the 
boundaries (diffusive boundaries) will decrease Tp as well 
and can potentially influence spin relaxation. Quantum 
mechanically, the channel narrowing leads to the quanti- 
zation of the electron transverse motion in the strip and 
absence of spin relaxation without intersubband scatter- 
ing. 

The first four possibilities are considered to some ex- 
tend in the scientific literature or are just evident conse- 
quences of the relaxation mechanisms. The fifth deserves 
a more thorough analysis. To check the effect of the pat- 
terning of the 2D electron gas into a strip of some large 
width L, we developed a simple Monte Carlo sumulator 
encompassing random scattering of the particles in the 
channel, reflection from the boundaries and spin rotation 



during free flight due to spin splitting of the conduction 
band. Now we describe in more detail our working model. 



III. MODEL 

As a model system we consider a strip of the 2D elec- 
tron gas. The third dimension (axis z) is quantizing and 
it is absolutely irrelevant to our consideration of the par- 
ticle movement in the real space, since the intersubband 
gap is larger than all other energy scales involved. 
i) We will assume that all particles have the same veloc- 
ity l^-l- 

ii) Scattering is considered to be elastic and isotropic in 
order to retain model simplicity; the former assumption 
preserves velocity modulus, the later one eliminates any 
correlations between directions of the particle velocity 
before and after the scattering event. 
Hi) We neglect electron-electron interactions and consider 
all particles to be independent. 

iv) An assumption that all scatterers are completely un- 
correlated leads to an exponential distribution of times 
between any two consecutive scattering events, their av- 
erage is called the momentum relaxation time, Tp. This 
time corresponds to the mean- free path Lp = \v\Tp. 
v) The problem spin Hamiltonian is given in the form of 
Eq. (0) and influences only spin coordinate. We ignore 
the reciprocal effect of the spin on the motion in the real 
space. 

vi) The width of the channel is large, several or even tens 
of mean-free paths so that it is permissible to consider 
classically the electron real-space movement. 
vii) At flrst we consider only reflecting boundary condi- 
tions at the borders of the 2D strip. Reflecting channel 
boundaries preserve longitudinal component of the par- 
ticle velocity and change sign of the normal component 
in collisions. Later, we will compare our results with dif- 
fusive boundaries. 



A. Types of experiments 

For simplicity, in all of our experiments particles will 
be injected into the system at some particular point A 
(input terminal) at time t — with spin 5'. As the 
particle experiences multiple scattering events a diffusive 
pattern of motion is formed with, for an isotropic system, 
a Gaussian distribution 



r(r) 



1 



■ exp 



Bp^/tfi 



that broadens as time increases: (r) '-^ ijp^t/Tp. 

Evidently, there are multiple possibilities for experi- 
mental setups. The definition of spin relaxation, obtained 
in these experiments will vary correspondingly. Now, we 
consider several important possibilities: 
i) The most informative type of experiment would be 



to obtain the average spin (S) as a function of time t 
at each point B (output terminal). Results of all other 
experimental configurations can be derived by partial in- 
tegration of this correlation function. 
ii) At time t average (5') is calculated for the whole 
ensemble independent of the real-space position of elec- 
trons. Optical experiments are considered likely to de- 
liver information of this type, because of the limited pos- 
sibilities to focus optical system and fundamental restric- 
tions. 

Hi) Particles, reaching output terminal are removed from 
the system immediately, (S) is measured as a function, 
i.e., of the interterminal distance. Individual particles 
can spend a substantial time in the system, depending 
on their trajectories. This type of experiment is the most 
probable variant in electric experiments where points A 
and B can be identified as real device gates. Made from 
ferromagnetic materials, gates can inject polarized elec- 
trons and sense the polarization of the drain flux, deliv- 
ering information about the average spin of carriers. 

Let us now show, that our result for spin relaxation 
does indeed depend on the deflnition of the experiment. 
As a simple example we consider a pure ID case. From 
the point of view of the first and third experiments there 
is no spin relaxation. There exists a systematic rotation 
of the spin proportional to the distance from the injec- 
tion point A. Indepent of the number of scattering events 
and individual trajectories, all particles reaching point B 
will have exactly the same spin orientation, but it would 
be different for different choices of point B. Thus, the 
transverse component of the spin Sy{x) = SyCOs{r]j:,px). 
For the second realization we readily conclude that the 
average spin 



{Sy) = / dX Sy{x)r{x) 



S'yCxp 



5°exp 



^VDpi^f 



1 2 r2 



that is an exponential decay of spin polarization in the 
spatially more and more broadening electron distribu- 
tion. 

Note that in the 2D case the same phenomenon of the 
systematic rotation of the electron spin takes place in ad- 
dition to the (2D~3D specific) DP spin relaxation. This 
systematic rotation on the angle 4> = yyop r x z for the 
particle real space transfer r is again independent on the 
detailes of individual trajectories. 



lation is performed for the D'yakonov-Perel' parameter 
VBpLp = 0.05 and different channel widths. For this par- 
ticular graph the trajectories of A^ = 5 x 10'^ electrons 
are traced that gives a standard deviation in the defini- 
tion of the spin polarization on the order of \/N ~ 10^^ 
(throughout our investigation A^ = 10"^ — 10^). 

A strong dependence of the polarization decay on the 
channel width is obtained. The decay is found to be 
approximately exponential, apart from the small initial 
interval. This is the region where the diffusive regime 
of electron motion and spin rotation is not yet estab- 
lished {^/t/Tp ~ 1). It overlaps or is followed by the 
region where the typical trajectory of a randomly scat- 
tered particle does not reaches channel boundaries yet. 
Here the relaxation should be essentially the same as in 
the unpatterned 2D gas. 

This generally exponential temporal behavior of spin 
polarization justifies an introduction of the spin relax- 
ation time, the dependence of this decay time on the 
channel width is presented in Fig. nb. In case of suf- 
ficiently wide channels so defined relaxation time ap- 
proaches the 2D limit of Tg^, for narrow channels ts 
scales as L~^. 

In Fig. nb we summarize results of our simulation 
of channels with a fixed width and different values of 
D'yakonov-Perel' term ?7dp- As the constant t^dp be- 
comes larger and larger, the system behavior graduately 
goes out of the regime of motional narrowing since ele- 
mentary rotations on the elementary free fiights are not 
small anymore. This leads to a quick spin relaxation as 
reflected by the reduction and saturation of rg. On the 
other side of this dependence is a steep increase in the 
relaxation time as 77DP decreases. We have found an in- 
termediate region, where Tg scales as a second power of 
the DP constant, that is foUowd by a forth-power depen- 
dence for very small values of 77DP • 

Thus, in the well-established regime of motional nar- 
rowing and in the limit of sufficiently narrow channels, 
we propose an asymptotic formula ts ^^ r/j^pL^^. 

To identify an effect of the particular choice of the 
boundary conditions we have simulated a channel with 
diffusive boundaries. The particle that reaches the chan- 
nel border is scattered back into the channel with equal 
probabilities for all directions of the new particle veloc- 
ity. The results of the calculation repeat closely the case 
of the reflective boudaries up to the very narrow channel 
with widths of only several mean free paths. No system- 
atic deviation of the spin relaxation time is observed for 
wider channels. 



IV. RESULTS OF SIMULATION 

Fig. |l|a shows the time dependence of the average spin 
polarization 2{Sy) in the channel. It is found that apart 
from the systematic rotation this dependence is essen- 
tially the same for all points inside channel that have a 
substantial electron occupation, so it is possible to cal- 
culate reliably a value of spin polarization. The calcu- 



V. REGIMES OF THE DP RELAXATION IN A 
CHANNEL 

Thus, we distinguish the following regimes of spin re- 
laxation in channels of finite width as we vary both in- 
dependent parameters yyop and L of the problem in con- 
sideration: 




5000 10000 15000 20000 

Time (in t) 



CO 10 r 



CD 

E 



X 

_C0 

0) 
C/D 10^ 



- ^ 


'7op;-,=0.05 


,-2 

s 

■ (b) 

1 


2D ■ ■ ■! 

^s ^ 

1 



10 

Width L (in /. ) 



100 




FIG. 1. Spin relaxation in a channel, (a) Time depen- 
dence of the spin polarization, calculated for different channel 
widths. DP constant rjop is fixed to be 0.05. The trajectories 
of 5 X 10^ particles are traced that defines a standard devia- 
tion ~ 10~^ for the calculated averages. Close-to-exponential 
decays of the polarization permit to define spin relaxation 
time, Ts; (b) rs as a function of the channel width L; (c) 
Spin relaxation time in dependence on the DP spin splitting 
constant t^dp at fixed channel width L = 40Lp. 
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FIG. 2. Different regimes of the spin relaxation on 

the plane (riuPjL) of model parameters: rjupLp > 1 
— elemenary rotations during free fiights are not small, 
Tp; rjupLp < 1, -qupL ^1 — 2D spin relaxation, 
— supression of spin relax- 

iHr T\-2 -4 r -2 

Ts (77DP-L) ~ Tp rij^pL ; 



TS 
Ts 



Tp{riupLp) ; rjupL < 1 
ation, quasi-lD regime, ts 
L<Lp 



L substitutes Lp, quantum mechanical quantiza- 
tion in the channel. 

Fh'st of all, for very large spin splitting {rj^pLp -^1) 
we violate a general condition for the motional narrow- 
ing regime for the DP spin relaxation. Each elementary 
rotation is not small and the information about the spin 
polarization is lost already after the first random scatter- 
ing event (see as well Fig. pi as a guide) . For this regime 
Ts ~ Tp and is the shortest of all regimes. 

When 77DP is small {rji^pLp < 1), we are back in the 
regime of motional narrowing and a well known equation 
Tg^ ~ Tp{'i]DpLp)^^ defines the time of spin relaxation in 
the 2D system. Now we narrow the strip of the 2D elec- 
tron gas. The behavior is unchanged until rj^pL '^ 1. For 
smaller channel widths {rfptpL < 1) DP spin relaxation is 
suppresed very effectively, with ts ~ t^^ {rjptp L)'"^ . For 
L < Lp the channel width L acts as a new mean-free 
path in the system, substituting Lp in equations. Actu- 
ally, this region does not satisfy our assumption about 
classical motion of particles in real space; the transverse 



motion is quantized and this system should be considered 
as a quantum wire with multiple subbands. The cross- 
over points define broad regions of mixed behavior that 
becomes more definitive as we move out of them. 

We conducted analysis presented here for spin relax- 
ation in the most convenient from our point of view 
units. Now, we will check that the chosen range of pa- 
rameters is well within the reasonable limits present in 
the contemporary heterostructures. For the asymmetric 
GaAs/AlGaAs quantum well with the electron concen- 
tration of 10^^ sm^^ in the channel the detected experi- 
mentally and confirmed by the calculation spin splitting 
in the conduction band is on the order of 0.2 — 0.3 meV 
atj-the Fermi energy (see recent paper by Wissinger et. 
alEj and references therein). This 2D electron concen- 
tration corresponds to the fcj? « 3 x 10^ cm~^. Sam- 
ples with Lp ^ 1/(7 = 3 X 10^^ cm are readily available 
in laboratories. Furthermore, this splitting corresponds 
to 77DP = 6 — 9 X 10^ cm^^. This set of parameters is 
on the border of the regimes of the motional narrowing 
and large elementary rotations. Reducing the spin split- 
ting (parameter ryop) (by, for example, reduction of the 
electron concentration in the 2D electron gas and, thus, 
structure asymmetry due to the internal electric field) 
will push parameters into well established regime of mo- 
tional narrowing, suitable for the experimental observa- 
tion of the described phenomena. Datta and DasQ give 
an estimate of Tr/ryDp = 7 x 10^^ cm for some particular 
InGaAs/InAlAs heterostructures. 



VI. SUMMARY 



In conclusion, we have investigated spin-dependent 
transport in semiconductor narrow 2D channels and ex- 
plored the possibility of suppressing spin relaxation. Our 
approach is based on a Monte Carlo transport model 
and incorporates information on conduction band elec- 
tron spins and spin rotation mechanisms. Specifically, 
an ensemble of electrons experiencing multiple scattering 
events is simulated numerically to analyze the decay of 
electron spin polarization in channels of finite width due 
to the D'yakonov-Perel' mechanism. We have been able 
to identify different regimes of the spin relaxation in the 
2D channels of finite width and obtain dependencies of 
the spin relaxation time on the width L and DP param- 
eter 77DP- The most attractive for the future spintronic 
applications is a regime of the suppresed spin relaxation 
with the relaxation time, ts, scaling as L~^. 
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